Quantum Optomechanics in the Bistable Regime 
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We study the simplest optomechanical system with a focus on the bistable regime. The covariance matrix 
formalism allows us to study both cooling and entanglement in a unified framework. We identify two key factors 
governing entanglement, namely the bistability parameter, i.e. the distance from the end of a stable branch in 
the bistable regime, and the effective detuning, and we describe the optimum regime where entanglement is 
greatest. We also show that in general entanglement is a non-monotonic function of optomechanical coupling. 
This is especially important in understanding the optomechanical entanglement of the second stable branch. 

PACS numbers: 



I. INTRODUCTION 

Observing quantum effects like superposition states or en- 
tanglement at the macroscopic level is a long standing goal. 
It is a widely held view that this should be possible, provided 
that enviromentally induced decoherence can be sufficiently 
suppressed. Note however that there are some theoretical pro- 
posals which would rule out the existence of quantum effects 
at the macroscopic level, see e.g. Ref. [1]. Proposals for the 
experimental observation of macroscopic quantum effects are 
often based on the principle of Schrodinger's cat, i.e. on cou- 
pling a microscopic quantum system to a macroscopic system 
in a controlled way, in order to create a macroscopic superpo- 
sition state 

One particularly promising approach in this context is the 
use of optomechanical systems. The most basic optomechan- 
ical system consists of a Fabry-Perot cavity with one mov- 
able end mirror. The position of this mirror is determined by 
the radiation pressure inside the cavity. Such systems were 
first studied in the context of high precision measurements 
and gravitational wave detection Jst] . It was suggested in Ref. 
||2il that the radiation pressure of a single photon in a high fi- 
nesse optical cavity could in principle create a macroscopic 
superposition of two spatially distinct locations of a move- 
able mirror . A potential implementation of this idea was 
proposed in Ref. |4]. It is very challenging experimentally 
to achieve sufficiently strong optomechanical coupling at the 
single-photon level, requiring a system that combines high op- 
tical and mechanical finesse, low mechanical resonance fre- 
quency and ultra-low temperature. 

One way to enhance the opto-mechanical interaction is to 
pump the cavity with a strong laser. Using this technique 
the strong coupling regime in optomechanical systems has re- 
cently been reached|6]. In the presence of a strong enough 
driving laser the field enhancement inside the high finesse op- 
tical cavity is large enough to trigger nonlinear behaviour of 
the system. Depending on the input power and the detuning 
of the driving laser with respect to the cavity resonance, op- 
tomechanical systems exhibit different types of nonlinear be- 
haviour. For strong enough inputpower, in the blue detuned 
regime one obtains multistabilityQ, instability (jsj] and chaotic 
motion |90. In the red detuned regime bistability lflo[[TTll oc- 



FIG. 1: Mechanical restoring force and radiation pressure force 
around a cavity resonance. The leftmost and rightmost intersection 
points are stable equilibrium positions, whereas the middle one is 
unstable. 



curs. Here we consider the red detuned regime. This is the 
appropriate regime for cooling the mechanical oscillator close 
to the ground state 112lfl3t], which is usually seen as a pre- 
requisite for observing quantum effects. We are particularly 
interested in the relationship between bistability and entan- 
glement. 

Optomechanical bistability can be understood intuitively as 
the result of a competition between the mechanical restor- 
ing force, which increases linearly when the mirror is moved 
from its equilibrium position, and the radiation pressure force, 
which has a maximum at the cavity resonance. For a suitable 
set of parameters as in Fig. 1 there are three intersection points 
between the two forces. The leftmost and rightmost intersec- 
tion points correspond to stable states, because the restoring 
force grows faster than the radiation pressure (for the right- 
most point the radiation pressure even decreases as the mir- 
ror is pushed outwards). In contrast, the middle intersection 
point is unstable because the radiation pressure force increases 
faster than the restoring force. 

It is known that the optomechanical interaction can squeeze 
the cavity mode, and this squeezing becomes maximal close 
to the bistable regime f\^. It has also been noted[ 15| that un- 
der certain conditions entanglement is maximized at the bista- 
bility threshold. This was interpreted as being due to the en- 
hanced optomechanical coupling strength in this region. Here 
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we analyze entanglement close to and in the bistable regime 
in detail. We show that a lot of insight can be gained by ana- 
lyzing the situation in terms of two key parameters, namely 
the effective detuning and the bistability parameter, which 
quantifies the distance from the end of each bistable branch. 
CooUng and entanglement can be studied in the same theo- 
retical framework based on the covariance matrix. We iden- 
tify the optimal regimes for both cooling and entanglement. 
We also show that, somewhat surprisingly, entanglement is in 
general a non-monotonic function of the optomechanical cou- 
pling strength. (Naively one might have expected it to always 
increase with optomechanical coupling strength .) 

The paper is organized as follows: Section II introduces the 
optomechanical system and describes the linearization of the 
equations of motion around the steady state. We also show 
how bistability arises in the red-detuned regime in this frame- 
work, introduce the bistability parameter, and derive the de- 
pendence of the photon and phonon number on this parameter, 
which leads us to a discussion of cooling. Section III discusses 
the optomechanical entanglement and its dependence on the 
bistability parameter and the effective detuning. This allows 
us to determine the optimum value for the detuning and the 
maximum achievable entanglement in our system. We discuss 
the role of the optomechanical coupling constant, show how 
entanglement varies on both stable branches in the bistable 
regime, and discuss its robustness under increasing tempera- 
ture. Section IV is a summary and conclusion. 



II. THE SYSTEM 

We consider a high Q Fabry-Perot cavity with decay rate 
K. The moveable mirror can move under the influence of ra- 
diation pressure and thermal noise. The moveable mirror is 
initially in equilibrium with a thermal bath at temperature T 
which results in the mechanical damping rate jm and the noise 
force ^{f). The system is driven by a laser with frequency ui 
and power P. The general Hamiltonian of such a system is 
derived in [16J . In the regime of parameters that we are inter- 
ested in, the general Hamiltonian simplifies to lfT6[ [TtIi 

H = hLO,a^a + ^^(p^+q^)-nGoa^aq+inE(a^e-''^'''-ae'"^'), 

(1) 

where ojc and a are frequency and annihilation operator of 
the cavity mode,respectively, w,„, q,p are frequency and di- 
mensionless position and momentum operator of the mir- 

ror,respectively, and Go = '^^e coupling constant 

and E - where P, oji are the input laser power and 

frequency respectively. The first two terms correspond to two 
free harmonic oscillators, the third term corresponds to the 
optomechanical coupling and the last term corresponds to the 
cavity being driven by the laser 

The equations of motion in the presence of damping and 
noise are 
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FIG. 2: Bistability of the intracavity power with respect to the in- 
put power. The solid and dotted lines correspond to the stable and 
unstable branches respectively. The inset shows the bistability pa- 
rameter T] for the two stable branches. The end of each stable branch 
corresponds to t; = 0. 

p = -a)„,q - y„,p + Gaa^a + ^(t) (3) 

a — -(k + ;Ao)a + iGoaq + E + VlKain (4) 

where Aq - oJc - w/,, «,„ is the vacuum input noise of the 
cavity, and ^(f) is the noise associated with the damping of 
the mechanical oscillator. The nonlinear Eqs. (3,4) can be 
linearized by expanding the operators around their steady state 
values Oi - Oi s + 60 i , where O, = a, q, p. 

From Eqs. (2,3,4), the steady state solutions are as = 

K+nAo^-Cocur " ^rlf^' " ^' ^^^'■^ a^^^p.are flie station- 
ary values for cavity amplitude, position and momentum of 
mechanical oscillator, respectively. Note that the last of these 
relations is a third order polynomial equation for as, which 
has three roots. The largest and the smallest roots are stable, 
and the middle one is unstable. Fig. 2 shows the hysteresis 
loop for the intracavity power Consider Pcav initially on the 
lower stable branch (I in Fig. 2, corresponding to the small- 
est root). As Pi„ increases past its value for the first turning 
point, Pcav switches to the upper stable branch (II in Fig. 2, 
coiTesponding to the largest root). For P,„ larger than its value 
at this switch point , Pcav is given by the upper branch. As 
Pi„ decreases below this value Pcav is still given by the upper 
branch until f ,„ decreases below its value for the upper turn- 
ing point. At this point Pcav switches back down to the lower 
branch. We have used the the set of parameters of ReffT?!, 
which is close to several optomechanical experiments [1 8-2 ij]. 
We consider Fabry-Perot cavity with length L = 1mm and fi- 
nesse 'F - 1.07 X 10"*, driven by a laser with A = 810nm and 
Ao = 2.62w,„. The mechanical oscillator frequency , damping 
rate and mass are lQMHz,lQQHz and 5ng respectively with 
environment temperature T = 400mA'. 

By introducing M^(f) = (6q(t),6p(t),X(t), Y(t)) and n^(f) = 
(0,^(0, ^Xi„(t), ^Yi„(t)) where X = ^^^^ and Y = 
^2^=2- and corresponding noises X,„and F,„, the linearized dy- 
namics of system can be written in a compact form 



q = cij„,p 



(2) 



u(t) - Au(t) + n(t) 



(5) 
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(6) 



and G — y/lGoOs , A = Aq - Goqs are the enhanced optome- 
chanical coupHng rate and effective detuning. 

Since the initial state of the system is Gaussian and the dy- 
namical equations are linear in the creation and annihilation 
operators both for cavity and mechanical mode the state of 
the system remains Gaussian at all times. A Gaussian state is 
fully characterized by its covariance matrix which is defined at 
any given time t by Vijit) = <"- W'o(')+'o(')", W> ^ mechanical 
and the optical input noise operators are fully characterized by 
their correlation function which in the Markovian approxima- 
tion are given by 



{ai„(t)al(t')) = 6(t - t'). 



= r„,(2«+ i)5(f-f'). 



(7) 



(8) 



where n - [ex/7( jy) - 1]"' is the mean thermal phonon num- 
ber and kB is Boltzmann's constant. From Eqs.(5,7,8) one ob- 
tains the equation of motion for the covariance matrix which 
is given by ll22ll 



V ^AV + VA' + D. 



(9) 



The steady state solution for the covariance matrix (V = 0) 
is reached if all the eigenvalues of the matrix A have negative 
real parts. In the red detuned regime of operation (A > 0), 
the Routh-Hurwitz criterion ll23ll gives the following stability 
condition 



G^A > 0. 



(10) 



In the following we use the dimensionless "bistability param- 
eter" defined as [24] 



^= 1- 



G^A 



W,„(/f2 + A2) 



(11) 



which is a positive number between zero and one according 
to Eq.(fTOli in the red detuned regime (A > 0). We have shown 
the bistability parameter in the inset of Fig. 2. As can be seen 
from Fig. 2, 77 decreases when approaching the bistable regime 
and becomes equal to zero at the end of each stable branch. 

Eq. (10) can be intuitively understood by ignoring retarda- 
tion eff'ects for the radiation pressure. Assuming that the op- 
tical field adiabatically follows the mechanical oscillator (i.e. 
setting 5a - 0), one has 5a - jj^^i-^^q + ^"I^Kain). The 
equation of motion for the mirror becomes 



5p = -(w,„ - 



G^A 



+ A2 

where ^7- = ^ + G ^Ki -^ 



)5q - y„,5p + 



(12) 



^.+,A ■ FromEq.(ll), wesee 

that the mechanical oscillator is stable if the first coefficient 



is positive. In this case the first term in Eq. (fT2T i corresponds 
to a harmonic restoring force, see also Fig. 1 and the associ- 
ated discussion. This implies the stability condition (10). The 
adiabatic approximation is equivalent to treating the response 
of the cavity field to the moving mirror as instantaneous. It is 
well known that the delayed nature of this response gives rise 
to cooling ['25'l, which is however not essential for the above 
argument. We feel that this argument helps the physical un- 
derstanding of the stability condition. However, let us empha- 
size that we will not make the adiabatic approximation in what 
follows. 

In the bistable regime the fluctuations around the steady 
state solution diverge as one approaches the end of each stable 
branch. To show this explicitly we solve the Eq.(9) for steady 
state, from which we can obtain the phonon and photon num- 
bers by using «,„ - and hg - Y^^lXazi^ xhe general 
solution is complicated and not very illuminating. Simple re- 
lations that show the dependence of the fluctuations on the 
stability parameter can be obtained by assuming a high me- 
chanical quality factor and low temperature environment,i.e. 

^ » 1 and ^ » 1 . We find 

r™ ny,,. 



Tlffi — 



8A?7fa;,„ 



8^ ■ 



(13) 



(14) 



From Eq.(13,14) it is clear that the phonon and photon num- 
bers diverge as rj approaches zero. In order to stay within the 
range of validity of the linearization approximation, we have 
made sure that So <sc I ffi P in all the results shown below. 

For rj ~ I from Eq.(13) the optimum value for the detuning 
which minimizes the phonon number is given by 



^opt — 



(15) 



Using the optimum value for detuning in Eq.(13) one finds 



"m = -( 1) 

z a>m 



(16) 



which is identical to Eq.(7) in 11261] . In the resolved sideband 
regime {a>,„ » k) one sees from Eq. (15) that ground state 
cooling can be achieved (n,„ - j^) ll26ll27ll . 



III. OPTO-MECHANICAL ENTANGLEMENT 



As shown in 112811. for bipartite Gaussian states the Peres- 
Horodecki criterion [TQ", '30) (positivity of the density matrix 
under partial transposition) is necessary and sufficient for sep- 
arability. In terms of the covariance matrix formalism this 
criterion is called logarithmic negativity which is defined as 



En - max{0,-ln(2v„in)}. 



(17) 
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FIG. 3: Optomechanical entanglement (a) and optomechanical cou- 
pling constant (b) as a function of bistability parameter and nor- 
malized effective detuning — , for a cavity decay rate k = 1.4w,„. 



where 



is the smallest symplectic eigenvalue of the 



partially transposed covariance matrix given by 



l-MlM^^ where E 



detA + detB - 2detC, and we repre- 
sent the covariance matrix in terms of 



A C 
B 



(18) 



Equipped with this measure we go on to study optomechan- 
ical entanglement. Fig. 3(a) shows the logarithmic negativity 
as a function of bistability parameter tj and effective detun- 
ing A. We note that for 77 ~ 1, which is required for ground 
state cooling, there is no optomechanical entanglement. From 
Figure 3(a) one can identify three different regimes depend- 
ing on the effective detuning. In the first regime, correspond- 
ing to A < Q.\(jL>m in the figure, there is no optomechan- 
ical entanglement. In the second regime, corresponding to 



0.1 w,„ < A < Q3(jJm in the figure, there is some optomechani- 
cal entanglement but the maximum value for entanglement is 
attained for values of the bistability parameter 77 somewhere 
between zero and one. This shows that for a fixed detuning, 
the maximum entanglement does not necessarily occur at the 
end of the bistable branch, cf. ifTsll . Finally in the third regime, 
A > Q3(jL>m in the figure, there is strong optomechanical entan- 
glement, and for each fixed value of detuning the maximum 
entanglement is in fact reached at the end of the branch, i.e. 
for 77 = 0. So in short we have Ej^ i < £^2 < £'a'3, where E^r, 
is logarithmic negativity in the /-th regime. 

Fig. 3(b) shows the corresponding optomechanical coupling 
in the different regimes. From Fig. 3(b) it is clear that optome- 
chanical coupling is monotonically decreasing function of ef- 
fective detuning i.e Gi > G2 > G3 where G, is the optome- 
chanical coupling constant in the /-th regime. So we see that in 
general entanglement is not a monotonically increasing func- 
tion of the coupling constant. These observations suggest that 
the key vaiiables that determine the entanglement behaviour 
are the effective detuning and the bistability parameter, not 
the optomechanical coupling constant. 

A more quantitative understanding of the different regimes 
for entanglement is possible by looking at the entanglement 
behaviour in the vicinity of 77 = 0. Assuming that — >> 1 



and » 1 one finds S 



a -H I and detV 



- , where 



b = 



A2 - + col 
16^ 



(A2 -I- /^2)(a2 SojI) 



c = 



2A\A^ + K^) + (A^ - k^)cjI 



128A4 



(A^ -I- /^0(4A^ + AA^K^ + ^A-'oji + ojD 
256A4w2 



(19) 



(20) 



(21) 



(22) 



From these equations it is possible to drive a simple form for 
logarithmic negativity. Close to the bistability region (77 << 1) 

we have - max{0, a + prj] where a - -ln(2 ^) and 

a _ (abd-b^c-d-) 
" ~ 2db^ 

It is worth noting that in contrast to the phonon and photon 
numbers, which diverge for 77 = 0, the logarithmic negativity 
has a finite limiting value given by a. While our linearization 
approximation is not justified for the point rj - itself, it 
does apply in its close vicinity, as the photon number drops 
precipitously as one moves away from the end point of the 
stable branch, cf. Eq. (14). 

Using our expression for entanglement close to bistability 
one can easily identify the three regimes shown in Fig. 3. The 
first regime corresponds to a,/3 < 0. The second regime cor- 
responds to ff < 0,j6 > or Q',j6 > and the third region 
corresponds to a > 0,y6 < 0. 
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FIG. 4: Plot of the optomechanical entanglement as a function of 
input power for both stable branches. The dot dashed (dashed) line 
corresponds to the end of the first (second) stable branch. The pa- 
rameters are the same as in Fie.2. 



Moreover as can be seen from Fig. 3 the maximum op- 
tomechanical entanglement is reached in the bistability region 
(for rj approaching 0) in the third regime. So the maximum 
achievable optomechanical entanglement is given by a. From 
Eqs .(19-22) the optimum value for effective detuning where 
entanglement takes its maximum value is 



^opt 




16(— )2-h81 



(23) 



Vox K - \A(jj,„ from Eq. (23) we obtain A^p, = 0.85(w,„. Com- 
paring this to Eq. (15) one sees that the optimum effective de- 
tuning values for cooling and entanglement are not the same. 
Even more importantly, the cooling performance is optimized 
for T] - \, whereas entanglement becomes maximal for rj - 0. 
Using the optimum value for detuning we obtain the follow- 
ing expression for maximum achievable entanglement in our 
system 



^N,max 



•In 



1 128^2 
5 V 8/f2 -H 45^2 



(24) 



Note that this takes its greatest possible value for k -0, giving 
£Mmax = -ln[3/5]=0.51. 

It is also interesting to look at the optomechanical entan- 
glement for the two stable branches and their behaviour in the 
bistable regime. Fig. 4 shows the logarithmic negativity as a 
function of input power for both stable branches. Varying the 
input power corresponds to varying rj, cf. Fig. 2. We note the 
persistence of entanglement in the second stable state in a very 
narrow window of parameter space. As can be seen in Fig.4 
the entanglement is maximum at the end of each branch, cor- 
responding to the third regime. The fast decreasing entangle- 
ment for the second branch is in agreement with the bistability 
parameter behaviour in Fig. 2. Inset shows the optomechan- 
ical coupling for different stable branches. One sees clearly 
that the coupling constant is not the decisive parameter for the 
amount of entanglement in our system, and in particular that 
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FIG. 5: Contour plot for bistability parameter (a) and entanglement 
(b) versus bare detuning Aq and input power P,,, in mW. The param- 
eters are the same as in Fig. 2 



the entanglement is a non-monotonic function of the coupling 
constant. 

Until now we studied the entanglement in terms of param- 
eters that are natural to use from a theoretical point of view. 
It is also interesting to look at entanglement in terms of pa- 
rameters that can be directly controlled experimentally. Fig. 
5 shows bistability parameter and entanglement as a function 
of bare detuning Aq and laser power P. Note that as we come 
close to the end of the branch for suitable detuning and large 
enough input power the entanglement increases. 

We have also studied the robustness of entanglement with 
respect to the temperature. The result is shown in Fig. 6 . One 
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FIG. 6: Plot of the logarithmic negativity versus the input power for 
different environment temperatures, T = OAK (solid line), T = 5 
(dot dashed), and T = 10 (dashed). The other parameters are the 
same as in Fig. 2. 

sees that for higher temperatures the entanglement survives 
only in the vicinity of the bistable region. 

Finally we note that in the recent experiment|61 the ratio 
of the input power to the critical power (the input power for 
which the bistability happen) is about 0.5. So the bistable 
regime should definitely be accessible experimentally. 

IV. CONCLUSION 

We have studied the simplest optomechanical system us- 
ing the covariance matrix formalism with a special emphasis 



on bistability. We recovered the standard results on optome- 
chanical cooling as a special case of our general expression 
for the phonon number However, our focus was on entangle- 
ment. We identified two key parameters, namely the effective 
detuning and the bistability parameter (i.e. the distance from 
the end of each stable branch in the bistable regime), and we 
showed that there are different regimes for entanglement as 
a function of these parameters. In particular we showed that 
maximum entanglement is achieved when the system is simul- 
taneously close to the red sideband (in terms of effective de- 
tuning) and close to the end of each stable branch (bistability 
parameter close to zero). We also showed that the dependence 
of entanglement on the optomechanical coupling is counter- 
intuitive, and that in the bistable regime the entanglement is 
particularly robust with respect to temperature increases . 

It would be very interesting to see experimental explo- 
rations of the phenomena described in the present paper How- 
ever, it should be noted that measuring the covariance matrix, 
which lies at the heart of our analysis, requires direct access 
to the position and momentum variables of the mirror, not just 
the quadratures of the light. It seems that this would require 
either an auxiliary measurement cavity, as proposed in Ref. 
UtI], or at least an additional laser beam. A detailed analysis 
of the resulting more complex dynamics is work for the future . 
In the present paper we focused on the entanglement charac- 
teristics of the basic system, which are already quite rich and 
intriguing. 
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